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ABSTRACT 

The connection between dark matter halos and galactic baryons is often not well-constrained nor well-resolved in 
cosmological hydrodynamical simulations. Thus, Halo Occupation Distribution (HOD) models that assign galaxies 
to halos based on halo mass are frequently used to interpret clustering observations, even though it is well-known 
that the assembly history of dark matter halos is related to their clustering. In this paper we use high-resolution 
hydrodynamical cosmological simulations to compare the halo and stellar mass growth of galaxies in a large-scale 
overdensity to those in a large-scale underdensity (on scales of about 20 Mpc). The simulation reproduces assembly 
bias, that halos have earlier formation times in overdense environments than in underdense regions. We find that the 
stellar mass to halo mass ratio is larger in overdense regions in central galaxies residing in halos with masses between 
® Mq. When we force the local density (within 2 Mpc) at z=0 to be the same for galaxies in the large-scale 
over- and underdensities, we find the same results. We posit that this difference can be explained by a combination 
of earlier formation times, more interactions at early times with neighbors, and more filaments feeding galaxies in 
overdense regions. This result puts the standard practice of assigning stellar mass to halos based only on their mass, 
rather than considering their larger environment, into question. 


1. INTRODUCTION 

In a A cold dark matter (ACDM) universe (Komatsu 
et al. 2011), the mass function and distribution of dark 
matter halos is well-understood. Dark matter only sim¬ 
ulations show that halos form in filaments that are con¬ 
nected at nodes containing (a) massive halo(s). This 
compares well with observations that show galaxies in fil¬ 
aments, sheets, clusters and superclusters (e.g. Springel 
et al. 2006). 

Connecting the stellar component of galaxies to dark 
matter halos can be difficult because galaxy growth de¬ 
pends on both internal and external factors. Gas accre¬ 
tion and mergers can drive galaxy star formation, while 
feedback from active galactic nuclei (AGN), supernovae 
and stellar winds can delay or stop growth. Many of these 
processes are difficult for cosmological hydrodynamical 
simulations to resolve, such as the radiative cooling of 
gas into molecular clouds, star formation, and energy in¬ 
fection into the interstellar medium from supernovae and 
AGN. 

Despite these complications, halo mass has been closely 
linked to a galaxy’s stellar mass and star formation. An¬ 
alytic models of galaxy formation through smooth gas 
accretion found that there is an upper mass threshold 
at which accreted gas is shocked to temperatures from 
which it can no longer cool within a Hubble time and file 
star formation, determining the mass at which galaxies 
must be red (Rees & Ostriker 1977; Silk 1977; Binney 
1977; White and Rees 1978). Recent numerical hydro¬ 
dynamical cosmological simulations by Keres et al. 2005 
(see also Ocvirk et al. 2008) have updated this model, 
showing that cold gas accretion can occur in halos with a 
range of masses, but the amount of cold gas in a galaxy is 
determined by its dark matter halo mass, Mhaio- There¬ 
fore, one would expect a galaxy’s Mhaio to be directly 
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linked to the SFR from cold gas accretion. This model is 
focused only on the central galaxies in halos, as these are 
the only galaxies for which we could expect a relation¬ 
ship between gas and dark matter halo mass based on 
the two-mode theory of gas accretion (Keres et al. 2005; 
Dekel & Birnboim 2006). A number of possible interac¬ 
tions can affect the gas, stellar, and dark matter mass 
of a satellite galaxy, as discussed in Boselli & Gavazzi 
(2006). 

Crain et al. (2009) use the GIMIC simulations- 
hydrodynamical “re”-simulations of the Millennium sim¬ 
ulation at (-2, -1, 0, -fl, -|-2 )ct of the mean density on the 
scale of ~20 Mpc-to determine whether halo mass is re¬ 
sponsible for the differences in the SFR density (SFRD) 
in different environments. They find that the SFR to 
halo mass ratio is the same in all environments, and con¬ 
clude that halo mass determines the rate at which galax¬ 
ies form stars. 

Because of the difficulty in directly modeling the bary- 
onic physics that drives gas cooling, star formation, and 
feedback, and using the theory that halo mass is the fun¬ 
damental parameter that determines other galaxy prop¬ 
erties, more analytic models have been used to connect 
galaxies to dark matter halos. In particular, standard 
Halo Occupation Distribution models (HODs) assume 
that halo mass is the fundamental parameter determining 
the stellar mass and internal processes of galaxies, and do 
not include any dependence on the larger environment in 
which the halo resides. These models are frequently used 
to make mock catalogues to interpret galaxy clustering 
measurements (e.g. Kauffmann et al. 1997; Jing et al. 
1998; Benson et al. 2000; Seljak 2000; Peacock & Smith 
2000; Ma & Fry 2000; Scoccimarro et al. 2001; Berlind 
& Weinberg 2002; Zheng et al. 2005; Tinker et al. 2008). 

Because the halo mass function depends on 
environment-more massive halos form in higher 
density regions of the universe, called halo bias (e.g. 
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Kaiser 1984)-these methods also reproduce the empir¬ 
ical finding that higher stellar mass galaxies reside in 
regions of higher galaxy density. For example, Abbas 
& Sheth (2006) compare galaxy clustering in the SDSS 
DR4 (M,r < -21) to measurements in mock catalogues 
created using HODs and to measurements from an 
analytic halo model calculation. The authors argue that 
the three samples agree well enough that correlations 
between galaxy properties can be entirely explained 
by the variation of the halo mass function in different 
environments (see also Skibba et al. 2006). Tinker et 
al. (2008) compare observed void statistics to those 
obtained using a standard HOD, and find that the sizes 
and emptiness of voids are in excellent agreement. 

However, more than the halo mass distribution has 
been shown to depend on the environment. Gao et al. 
(2005) use the Millennium Simulation to find that the 
ACDM paradigm predicts that the clustering of dark 
matter halos depends not only on their mass but also 
on their formation time (the time at which the halo 
mass of the main progenitor has reached half of the final 
halo mass), an effect often called assembly bias. Specifi¬ 
cally, low-mass halos (Mhaio < 6.15x10^^ h~^ Mq) that 
assemble early are much more strongly clustered than 
those that assemble late (see also Harker et al. 2006; 
Gao & White 2007; Wechsler et al. 2006; Wetzel et 
al. 2007). Glustering, measured using the 2-point au¬ 
tocorrelation function, increases with increasingly early 
formation times, and the signal is stronger as the length 
scale increases. Fakhouri & Ma (2009; 2010) find that 
the formation time of dark matter halos in the Millen¬ 
nium Simulation is earlier with increasing environmental 
density, measured using the dark matter mass within 7 
Mpc of a galaxy. They further find that mass growth 
proceeds differently in different environments, with more 
mergers in high-density regions and more diffuse accre¬ 
tion in low-density regions. As Gao et al. (2005) point 
out, galaxy properties may depend on the assembly his¬ 
tory of halos, and therefore models that ignore the age 
dependence of clustering do so at their peril. 

Berlind et al. (2006) examine SDSS groups and find 
that central galaxy color is correlated with clustering, 
but only for the most massive (>10^^ M 0 ) galaxy groups. 
Groups with less red central galaxies cluster more than 
groups with redder central galaxies. They conclude that 
massive halos that formed earlier contain redder galaxies 
than more recently-formed halos. However, Yang et al. 
(2006) use the emission- and absorption-line strength in a 
galaxy’s spectrum to study the SFR of the central galax¬ 
ies of SDSS groups. They find that group halos at all 
masses with central galaxies with lower SFRs are more 
clustered than groups with highly star-forming central 
galaxies. Although Berlind et al. (2006) point out that 
emission- and absorption-line strength in a luminous red 
galaxy does not necessarily reflect its g — r color, it is not 
clear why these two groups get these seemingly opposing 
results. 

Much work has been performed studying assembly bias 
using SAMs or HODs. For example, Gonzalez & Padilla 
(2009) use a SAM to study the effect of the environment 
on galaxy SFRs. They find that most of the differences 
in the SFR between galaxies in different large-scale envi¬ 
ronments can be explained by the environmental depen¬ 


dence of the halo mass function. Further, they conclude 
that assembly bias is the most likely candidate driving 
the small differences in the populations beyond those ex¬ 
plained by the halo mass function. 

Groton et al. (2007) come to a different conclusion us¬ 
ing a SAM built into the Millennium Simulation. They 
compare the two-point autocorrelation functions for the 
galaxies formed in their SAMs to a population that has 
been shuffled randomly into halos of the same mass 
(maintaining the relative separation between the central 
and its satellites). Shuffling reduces the clustering of the 
entire sample. They then change their shuffling scheme 
to also account for halo formation time (t 5 o%) or concen¬ 
tration in addition to halo mass, but find that neither of 
these halo properties can account for the clustering dif¬ 
ferences seen in the SAM population versus the shuffled 
population. They conclude that another unknown aspect 
of halo assembly or environment must be causing a large 
fraction of the bias. 

Zentner et al. (2014) find that assembly bias affects 
HOD modeling by comparing the HOD fits to mock cat¬ 
alogues with and without assembly bias. They find that 
reasonable levels of assembly bias in the population can 
lead to systematic errors in the galaxy-halo connection 
inferred using standard HOD models, and conclude that 
incorporating assembly bias effects into future HODs 
should be a priority. 

Jung, Yee & Li (2014) use a SAM to study whether 
the stellar mass of galaxies with the same halo mass is 
affected in regions of different large-scale (7 h~^ Mpc) 
density. They find a small difference, specifically that 
low-mass (<10^^ Mq) halos have slightly higher stellar 
masses in high-density environments. 

In this paper we compare halo mass growth, stellar 
mass assembly and the star formation histories of galax¬ 
ies from z=6 to z=0 in a large-scale (^20 Mpc) over¬ 
density and underdensity cosmologically simulated using 
the adaptive mesh rehnement (AMR) Eulerian hydrody- 
namical code Enzo (Bryan et al. 2014). At their respec¬ 
tive volumes, they represent -|-1.8cr and -l.Ocr fluctua¬ 
tions. These simulations provide us with an opportunity 
to compare different large-scale environments within the 
same fully hydrodynamical cosmological simulation per¬ 
formed at sub-kpc resolution. It is important to compare 
across different environments within the same simula¬ 
tion, as all physical processes are modeled in the same 
way. 

Our goal is to determine whether the stellar mass of 
galaxies is universally related to the halo mass of galax¬ 
ies, or if it also depends on environment. After a brief 
description of our simulations (Section [^l, we discuss our 
galaxy selection technique and our method for determin¬ 
ing quantities for each galaxy in Section l2.ll In our re¬ 
sults we compare the formation histories of galaxies in 
different environments, binned either by their final halo 
mass (Section [2T]) or stellar mass ('Section [3.21) . In Sec¬ 
tion |4] we discuss several possible causes for the higher 
stellar masses found in galaxy halos in the overdense en¬ 
vironment, ending with our comprehensive interpretation 
in Section lTTl We then compare our results to other the¬ 
oretical work on this issue (Section[5]). Finally, in Section 
E we summarize our conclusions and discuss the implica¬ 
tions of our results. 
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2. METHOD 

For the details of our simulations, we refer the reader 
to Cen (2012), although for completeness we reiterate 
the main points here. We perform cosmological simula¬ 
tions with the AMR Eulerian hydrodynamical code Enzo 
(Bryan 1999; Bryan et al. 2014). We use cosmologi¬ 
cal parameters consistent with the WMAP7-normalized 
LCDM model (Komatsu et al. 2011): JIm = 0.28, fit, 
= 0.046, flA = 0.72, CT8 = 0.82, Ho = 100 h km s'^ 
Mpc“^ = 70 km s“^ Mpc“^, and n = 0.96. We first ran 
a low resolution simulation with a periodic box of 120 h~^ 
Mpc on a side, and identified two regions: an overden¬ 
sity centered on a cluster and an underdensity centered 
on a void at z = 0. We then resimulated each of the 
two regions separately with high resolution, but embed¬ 
ded within the outer 120 h~^ Mpc box to properly take 
into account large-scale tidal field effects and appropri¬ 
ate fluxes of matter, energy and momentum across the 
boundaries of the refined region. 

The overdense refined region, or C box, is 21 x 24 x 
20 h~^ Mpc^. The central cluster has an M 200 of ~2 x 
10^"^ Mq with a virial radius (r 2 oo) of 1.3 h~^ Mpc. The 
underdense refined region, or V box, is somewhat larger, 
at 31 X 31 X 35 h~^ Mpc^. At their respective volumes, 
they represent -|-1.8cr and -I.Oct fluctuations where a is 
the density variance on the volume of the C and V boxes. 
Although these are large-scale over- and underdense en¬ 
vironments, these high-resolution boxes are much larger 
than the cluster or the void at their centers. Thus, there 
are galaxies at a range of local densities in both boxes, 
and there is substantial overlap of local densities between 
the two volumes (Tonnesen & Cen 2012). 

In both refined boxes, the minimum cell size is 0.46 
h~^ kpc, using 11 refinement levels at z = 0. The initial 
conditions for the refined regions have a mean interpar¬ 
ticle separation of 117 h~^ kpc comoving, and a dark 
matter particle mass of 1.07 x 10® h~^ Mq. 

The simulations include a metagalactic UV back¬ 
ground (Haardt & Madau 1996), a model for shielding 
of UV radiation by neutral hydrogen (Cen et al. 2005), 
and metallicity-dependent radiative cooling (Cen et al. 
1995). The fraction and density of neutral hydrogen is 
directly computed within the simulations. Star particles 
are created in gas cells that satisfy a set of criteria for 
star formation proposed by Cen & Ostriker (1992), and 
reiterated with regards to this simulation in Cen (2012). 
Briefly, A star particle is created if the gas in a cell at 
any time meets the following three conditions simulta¬ 
neously: (1) contracting flow, (2) cooling time less than 
dynamical time, and (3) Jeans unstable. A star particle 
of mass m* = is created (the same amount 

is removed from the gas mass in the cell), where At is the 
time step, U = max(tdj^„,10’^ yr), tdyn = \//{32Gptot) 
is the dynamical time of the cell, nigas is the baryonic 
gas mass in the cell, and c* = 0.03 is the star formation 
efficiency. Each star particle has a mass of ^10® Mq, 
which is similar to the mass of a coeval globular cluster. 
Once formed, the stellar particle loses mass through gas 
recycling from Type II supernovae feedback, and about 
30% of the stellar particle mass is returned to the ISM 
within a time step. Supernovae feedback is implemented 
as described in Cen (2012): feedback energy and ejected 
metal-enriched mass are distributed into 27 local gas cells 


centered at the star particle in question, weighted by 
the specific volume of each cell, which is to mimic the 
physical process of supernova blastwave propagation that 
tends to channel energy, momentum, and mass into the 
least dense regions (with the least resistance and cool¬ 
ing). We allow the whole feedback process to be hy- 
drodynamically coupled to surroundings and subject to 
relevant physical processes, such as cooling and heating, 
as in nature. The simulation used in this paper has com¬ 
pared several galaxy properties that depend critically on 
the feedback method to observations and found strong 
agreement (Cen 20IIa-c; Cen 2012). 

We do not include a prescription for AGN feedback in 
this simulation, and as a result, our simulation overpro¬ 
duces luminous galaxies at the centers of groups and clus¬ 
ters of galaxies. This is discussed in detail in Cen (2011c), 
who shows that the luminosity function of the simulated 
galaxies agrees well with observations at z=0 except at 
the high-luminosity end. When Cen (2011c) adds an 
AGN feedback correction in the post-simulation analysis 
that most strongly affects halos with masses greater than 
10^® Mq or galaxies with stellar masses above 4x10®^ 
M 0 , the simulated luminosity function also agrees with 
observations at the high luminosity end. We do not in¬ 
clude any post-simulation AGN correction because of un¬ 
certainties in its implementation, particularly across a 
large redshift range (see discussion in Tonnesen & Cen 
2014), and as all of the galaxies we examine in this paper 
have halo masses less than 10^® Mq, we would expect the 
correction to be minor in any case. 

2.1. Galaxies 

We use HOP (Eisenstein & Hut 1998) to identify galax¬ 
ies using the stellar particles. HOP uses a two-step proce¬ 
dure to identify individual galaxies. Eirst, the algorithm 
assigns a density to each star particle based on the distri¬ 
bution of the surrounding particles and then hops from 
a particle to its densest nearby neighbor until a maxi¬ 
mum is reached. All particles (with densities above a 
minimum threshold, 6outer) that reach the same maxi¬ 
mum are identified as one coherent group. The densest 
star particle is considered the center of the stellar group, 
or galaxy. In the second step, groups are combined if 
the density at the saddle point which connects them is 
greater than 5saddle- The minimum number of star par¬ 
ticles in a group that HOP will return as a galaxy is 20, 
but this lower mass limit of 2x10^ Mq is well below the 
lowest mass of the galaxies followed in this paper at z=0, 
which is well above 10® Mq. We use HOP because of its 
physical basis, although we expect similar results would 
be found using a friends-of-friends halo finder. HOP has 
been tested and is robust using reasonable ranges of val¬ 
ues (e.g. Tonnesen, Bryan, & van Gorkom 2007). 

In this paper we measure the stellar mass (M*), dark 
matter halo mass (Mhaio), and SFR of galaxies. Stellar 
mass is determined by adding the mass of each star par¬ 
ticle identified by HOP to belong to a galaxy. The dark 
matter halo mass is calculated by summing the mass of 
all the dark matter particles out to r 2 oo (the radius within 
which the average density of the dark matter halo is 200 
times the critical density). 

We also categorize our galaxies into centrals and satel¬ 
lites. Central galaxies are those galaxies that are not 
within r 2 oo of any more massive galaxy, using the stel- 
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lar mass as identified with HOP. Satellite galaxies are 
defined as within r 2 oo of a more massive galaxy. Our 
qualitative results do not depend on our specific choice 
of radius, and are qualitatively the same if we use a ra¬ 
dius of 2 r 2 oo- In this paper, because we are studying the 
connection between halo mass and stellar mass, we only 
focus on central galaxies. 

Lackner et al. (2012) describe in detail the merger trees 
that we will use in this paper, updated to include addi¬ 
tional redshift outputs. Briefly, the method tracks the 
stellar particles from one galaxy to another in sequential 
redshift slices. If most of the stars in galaxy B at red- 
shift 0.05 are in galaxy A at redshift 0, then we consider 
galaxy B to be the parent of galaxy A. If several galaxies 
contribute stars to galaxy A, the galaxy that contributes 
the most stellar mass to galaxy A is the main progen¬ 
itor. A galaxy can be tracked to the minimum stellar 
mass identified by HOP, ~2xl0^ Mq, and this sets the 
minimum initial mass possible for galaxies we track. We 
note that not all galaxies are initially identified with such 
a small set of particles. 

In this paper we follow the stellar and dark matter 
mass histories of galaxies, so select only galaxies that 
have been tracked through the merger trees from a mini¬ 
mum redshift. As we are examining whether galaxy stel¬ 
lar mass grows differently in low- versus high-density en¬ 
vironments, we focus on galaxies that can be tracked to 
at least a redshift of I. We compared these results to a 
sample that could be tracked to at least a redshift of 3 to 
include the peak of the SFRD (Hopkins & Beacom 2006; 
Karim et al. 2011; Seymour et al. 2008; Brinchmann et 
al. 2004; Reddy & Steidel 2009; Bouwens et al. 2007), 
and found no qualitative differences in our results. 

In fact, we also find that if we compare the entire set of 
HOP-identified galaxies in each environment (as in Sec¬ 
tion [5]), we find the same results as when we focus on our 
tracked samples of galaxies. This is an important ver¬ 
ification exercise because demanding that a galaxy can 
be tracked to z=l is essentially demanding that a galaxy 
has a stellar mass of at least 2x10^ Mq by that redshift. 
Particularly in our lowest mass bin, this will tend to se¬ 
lect galaxies with higher z=0 stellar masses. Because we 
find the same results when we include all galaxies at in¬ 
dividual outputs, we know that this selection criterium 
does not effect our results. However, as one of the goals 
of this paper is to consider the growth of galaxies over 
time, we focus on the tracked sample. 

3. RESULTS 

3.1. Matching Halo Mass 

In order to compare the properties of galaxies over 
time, we have two galaxy samples: those in the large- 
scale underdensity (V box) and those in the large-scale 
overdensity (C box), both of which must be tracked to at 
least z=l. Recall that in this paper we only consider cen¬ 
tral galaxies. While this includes most galaxies in the V 
box (more than 90% of galaxies at z=0), at lower masses 
in the C box there are a large fraction of satellite galax¬ 
ies (~35% at 2=0), so the mass growth of the galaxies 
we follow is not necessarily representative of the C box 
population as a whole. In order to present our results the 
most clearly, we have binned both samples by halo mass 
at 2=0. We have verified that the specific mass ranges 


do not affect our results: we have selected our samples in 
several different ways and found that our results are ro¬ 
bust, even when we use C box samples that have slightly 
lower halo masses than their V box counterparts. We 
consider four mass bins, with Mhaio/^Q between 10^^- 
1011.5, 10“-9-10i2-3, and 

In Figures [T][5] we show several properties of our galax¬ 
ies over our tracked redshift range. In all figures, the solid 
lines are the median value of the sample, and the shading 
denotes the middle range of values (25-75 percentiles). 
The V box galaxies (underdense environment) are shown 
in blue and the C box galaxies (overdense environment) 
are shown in red. The vertical lines are the median red¬ 
shift at which the halo has grown to half of its final mass 
(formation time). Until our minimum tracking redshift, 
2 = 1 , new galaxies can be added to the set, which may 
decrease the median halo or stellar mass values. This can 
result in the sometimes dramatic differences in median 
values as galaxies are added to the sample, particularly 
in the lowest mass bin. These jumps are not important 
to our results, which as we have discussed are insensitive 
to our sample selection, including the redshift to which 
we track galaxies. 

In order to speed post-processing of our galaxy param¬ 
eters, we use coarse radial bins to determine r 2 oo, and do 
not calculate the average overdensity continuously. This 
means that the halo mass can vary, as seen particularly 
at late times in the underdense box, because a small dif¬ 
ference in the radial density profile can result in a large 
variation in our r 2 oo due to the bin size, and therefore 
Mhaio may vary. These variations also do not affect our 
results. 

3.1.1. Halo Mass 

First we consider the halo mass of our galaxy samples 
in Figure[TJ We can compare the formation time (or halo 
age) of our tracked galaxies in the overdense and under¬ 
dense environments. As in Gao et al. (2005), we define 
the age of a halo using the redshift at which the mass of 
the main progenitor halo is 50% of the final halo mass. 
Our outputs are in redshift steps of 0.05 and 0.1 in the 
C and V box, respectively, which sets our uncertainty 
on the halo age of any individual galaxy. If we do not 
start tracking a galaxy until after it has reached half of 
its final mass, we set the formation time to the earliest 
time tracked. We compare our median formation time 
including these estimated values to a median formation 
time calculated using only the formation times for galax¬ 
ies that we begin tracking when they are less than 50% of 
their final halo mass, and find that for the three highest 
bins the change is minuscule. The lowest bin, that tends 
to begin tracking galaxies at later times, shifts to earlier 
median formation times when only the smaller sample is 
used. In the lowest mass bin, from Mq, the 

formation times go from being similar in the two envi¬ 
ronments to being at higher redshifts in the overdense 
sample. 

If we compare across different masses in the same en¬ 
vironment, we find that in our underdense samples, as 
in Behroozi et al. (2013), the halo age decreases as halo 
mass increases. This is evident when comparing the high¬ 
est mass bin to any of the lower mass bins. In contrast, 
the overdense galaxy samples tend to have slightly de¬ 
creasing formation times with increasing mass until the 
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Fig. 1.— The median halo mass of tracked galaxies from a pos¬ 
sible maximum redshift of six. Galaxies are binned according to 
halo mass at 2=0. The solid lines are the median value of the sam¬ 
ple of tracked galaxies, and the shading denotes the middle range 
of values (25-75 percentiles). The V box galaxies (underdense en¬ 
vironment) are shown in blue and the C box galaxies (overdense 
environment) are shown in red. The vertical lines are the median 
redshift at which the halo has grown to half of its final mass. Our 
galaxy sample selection has selected gala xies with similar 2=0 halo 
masse s in the V and C boxes. See Section rS.ll for details and Section 
l3.1.1l for discussion. 

most massive halo bin ® M 0 ), which has the 

oldest formation time. It is possible that our selection 
criterion that a galaxy must remain a central galaxy to 
z=Q eliminates the older low-mass galaxies, which will ef¬ 
fect the median formation time in the overdense environ¬ 
ment more strongly because a higher fraction of galaxies 
are satellites. 

Comparing between galaxies of the same halo mass in 
underdense versus overdense environments, we find that 
in general halos in the overdense region have older or 
similar formation times than galaxies in the underdense 
environment. This tendency agrees with the halo bias 


Fig. 2.— The median stellar mass of tracked galaxies binned 
by halo mass at redshift 0. Lines are as in Figure ^ The C box 
galaxies in the overdense environment (red) tend to have higher 
final stellar masses than the V box galaxies in the underdense en¬ 
vironment (blue), and the crossover redshift from lower to higher 
stellar masses increases with increasing redshift. See Section l3.1.2l 
for discussion. 

findings of Gao et al. (2005), that more clustered halos 
tend to collapse first. As we expect, the biggest differ¬ 
ence is in the highest mass bin (lO^^'^-lO^^ ® M©), where 
we are the least likely to exclude halos from our sample 
because they become satellites. 

We find that, in agreement with the similar halo ages, 
the halo growth curves are very similar in the lower mass 
bins (lO^^-lO^^'^ Mq and lO^^'^-lO^^ ® Mq), particularly 
after about 2=2.5. However, in the highest Mhaio bin, 
the C box galaxies’ halo growth flattens from about 2=1, 
while the galaxies in the V box continue to grow more 
massive. 

In summary, in agreement with previous N-body sim¬ 
ulations (e.g. Gao et al 2005 and see our Introduction), 
we find evidence of assembly bias. Also, in the highest 
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Fig. 3.— The median stellar mass to halo mass ratio (SMHM) 
of tracked galaxies binned by halo mass at redshift 0. Lines are 
as in Figure^ The galaxies in the overdensity (C box, red lines) 
have higher SMHM ratios than those in the underdensity (V box, 
blue lines). This result holds e ven com paring between galaxies in 
different mass bins. See Section 13.1.21 for discussion. 

mass bin ® Mq), we see a hint that at low 

redshift (2:<1) the dark matter mass growth of massive 
halos flattens in galaxies in the overdense environments 
in comparison to halos in the underdense environment. 

3.1.2. Stellar Mass 

In Figure[2]we focus on the stellar mass of our galaxies. 
In the overdense (C box) environment, the stellar mass 
growth from z=l to z=0 is the steepest in the lowest halo 
mass bin (lO^^-IO^^'® Mq), indicating a trend of downsiz¬ 
ing in stellar mass growth with decreasing redshift. In all 
mass bins in the underdense environment and in the up¬ 
per three mass bins in the overdense sample the stellar 
mass growth from z<l is similar. However, the galax¬ 
ies in the underdense environment grow more than the 
galaxies in the overdense environment in all but the low- 


Fig. 4.— The median SFR of tracked ralaxies binned by halo 
mass at redshift 0. Lines are as in Figure |l] The SFRs of galaxies 
in the overdensity (C box, red lines) tend to be lower than those 
of galaxies in the underdensity (V box, blue lines ) at z'>A, but 
become higher at lower redshifts. See Section 13.1.31 for discussion. 

est mass bin, indicating that stellar growth shifts from 
overdense to underdense regions over time. 

In all of the halo mass bins, the median stellar mass 
of galaxies in the overdense box is higher than that of 
galaxies in the underdense box by z=0. In all but the 
highest mass bin ® Mq), the median stellar 

mass in the C box galaxies begins lower than the median 
stellar mass for the V box galaxies. The redshift at which 
the C box stellar mass crosses the V box stellar mass 
tends to decrease with decreasing mass-from the highest 
mass bin always having higher stellar mass galaxies in 
the C box to a crossover at z~I in the lowest mass bin. 
However, the stellar mass growth at very low redshift 
(z<I) is generally flatter in the C box galaxies than in the 
V box galaxies, indicating less recent star formation, in 
agreement with observations (e.g. Szomoru et al. 1996). 

The stellar mass to halo mass (SMHM) ratio, shown 
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Fig. 5. — The median sSFR of tracked galaxies binned by halo 
mass at redshift 0. Lines are as in Figure^ Galaxies follow the 
sSFR-stellar mass relation. Comparing galaxies with the same z=0 
halo mass, we find that at early times the sSFR of galaxies in the 
overdensity (C box, red lines) is higher with respect to the sSFR of 
galaxie s in th e underdensity (V box, blue lines) than at z=0. See 
Section l3.1.3l for discussion. 

in Figure [3l also reflects the higher stellar mass in halos 
in the overdense environment. In general, even when 
comparing across different halo mass bins, we see that 
the SMHM ratio is higher in the overdense environment, 
showing that the halo mass distribution does not effect 
this result. In other words, by z=0, the SMHM ratio 
of overdense galaxies in any mass bin is higher than the 
SMHM ratio of underdense galaxies in any halo mass bin. 

Comparing within each environment, we find that the 
SMHM ratio at z=0 is the highest for galaxies with the 
highest halo masses, unlike previous results (e.g. Moster 
et al. 2012, Guo et al. 2010, Behroozi et al. 2013 and 
references therein) that the peak SMHM ratio will be 
found for galaxies with halo masses of about 10^^ Mq. 
These previous results use several methods to match re¬ 


sults to existing data, and we do not call their overarch¬ 
ing SMHM results, that do not consider different large- 
scale environments, into question. Indeed, our SMHM 
ratio tends to be too high. This may indicate that our 
feedback mechanisms are inefficient, particularly at the 
lowest and highest masses, which we will discuss in de¬ 
tail below in Section Mi We conclude that the low¬ 
est mass and highest mass halos have the most cold gas 
that is able to form stars, either by never shock-heating 
accreting gas to high temperature or by having a high 
enough central density to cool gas and form stars. As 
our feedback prescription is the same across all environ¬ 
ments, the differences between galaxies in the overdense 
and underdense regions should be real and not due to 
our particular numerical methodology, although we can¬ 
not eliminate the possibility that our specific feedback 
scheme could smooth out or exacerbate any differences 
between galaxies in different environments. 

We conclude that the SMHM ratio differs in these 
differing large-scale environments, with central galaxies 
in overdense environments having higher SMHM ratios 
than those of central galaxies in underdense environ¬ 
ments. 

3.1.3. Star Formation Rate 

We next focus on the SFRs of the galaxies in Figure IH 
For all galaxy samples but the highest mass galaxies in 
the underdense environment, the SFR initially increases 
to an early peak, then decreases to z=0. This follows our 
expectations, because as shown in Cen (2011c) our sim¬ 
ulated SFR density history agrees well with observations 
that show a peak between z=2-3 followed by a declining 
SFRD. 

When we compare galaxies of the same mass in dif¬ 
ferent environments, we find that for all but the highest 
mass galaxies, the SFR in galaxies in the overdense en¬ 
vironment begins lower than the SFR in galaxies in the 
underdense environment. The SFR of the C box galaxies 
increases with respect to the SFR of the V box galaxies 
as the universe ages, eventually overtaking the SFR of V 
box galaxies. The redshift at which the C box SFR be¬ 
comes higher than the V box SFR tends to decrease for 
decreasing halo mass. In our highest mass bin (10^^'^- 
10^2 ® M 0 ), the C box SFR begins higher than the V box 
SFR, but decreases quickly after z=l. 

In Figure [5] we see that in either environment at z=0 
our galaxies follow the sSFR-stellar mass relation, with 
lower mass galaxies having higher sSFRs. We also see 
that at early times the sSFR of galaxies in the overdense 
environment tends to be somewhat higher than the sSFR 
of galaxies in the underdense environment. At late times, 
the z=0 sSFR is either similar in the two environments, 
or higher in the V box than in the C box, probably re¬ 
flecting the sSFR-M* relation. This may also agree with 
the Szomoru et al. (1996) observational result that galax¬ 
ies in the Bootes Void tend to lie above the Tully-Fisher 
relation using the B-band magnitude. 

The comparison between the star formation histories 
of these galaxies align with the comparative stellar mass 
growth of these galaxies. 

3.2. Matching Stellar Mass 

In order to compare the properties of galaxies using a 
more easily observed quantity, we verify our results when 



















we bin galaxies by their stellar mass at z=0. We find 
that, as when we bin by halo mass, galaxies in the over- 
dense region tend to have higher SMHM ratios. Unlike 
when we bin by halo mass, we find that the SFRs and sS- 
FRs in the two environments are more similar by z=0, as 
we might expect given that the final stellar mass is quite 
similar. However, the z=0 SFR and sSFR of our mas¬ 
sive galaxies in the underdense region tends to be slightly 
higher than in the overdense region. Thus, depending on 
the stellar masses of the galaxies we compare, the sSFR is 
the same across different environments, or galaxies of the 
same mass are bluer in underdense environments than in 
overdense environments (Yang et al. 2006; Rojas et al. 
2004; Grogin & Geller 1999,2000; Szomoru et al. 1996). 

4. DISCUSSION 

We will now discuss several possible causes of this dif¬ 
ference in the SMHM ratio in galaxies in different envi¬ 
ronments. Finally, we will propose our synthesized pic¬ 
ture of how higher-density large-scale environments could 
produce galaxies with high SMHM ratios. 

4.1. Formation Time 

We first discuss whether the formation time, defined 
as in Gao et al. (2005) and marked with vertical lines in 
Figures [T][^ is related to the increased SMHM ratios in 
the overdense environment. We see a clear trend that the 
C box galaxy populations tend to have earlier formation 
times and higher SMHM ratios. Although the sample 
used in this paper of overdense galaxies is slightly more 
massive than the galaxies in the underdense region, we 
have verified that this result does not depend on the mass 
distribution of the two samples by selecting an overdense 
galaxy sample that has slightly lower masses than the 
underdense sample. 

However, we hnd that more subtle differences in for¬ 
mation time are not clearly related to the final SMHM 
ratios. Examining Figure [31 we see that the two largest 
differences in the SMHM ratios are in the halo mass bins 
from - 10^^ ® and 10^^ ® - 10^^'^. These mass bins 

have the smallest differences in the formation times be¬ 
tween the overdense and underdense populations. The 
largest difference in formation times is in the highest 
mass bin ® Mq), and the difference between 

the hnal SMHM ratios in the two environments is among 
the smallest. 

To summarize, we find that galaxies in the over- 
dense environment have both earlier formation times and 
higher SMHM ratios than galaxies in the underdense en¬ 
vironment. However, an earlier formation time does not 
necessarily lead to a higher SMHM ratio, as seen by 
comparing the galaxies in the underdense environment 
in Figure [31 and, as we discussed above, a larger differ¬ 
ence in formation times does not force a larger difference 
in SMHM ratio. 

4.2. Tidal Stripping of DM Halos? 

Several recent works have argued that assembly bias, 
the earlier formation times of halos in higher density en¬ 
vironments, can be explained by the increased tidal in¬ 
teractions with neighbors in higher density environments. 
N-body simulations have shown that at late times, the 
tidal field from massive neighbors can halt halo growth 


and even in some cases reduce halo mass (Avila-Reese 
et al. 2005; Maulbetsch et al. 2007; Wang et al. 2006; 
Diemand et al. 2007; Hahn et al. 2009; McBride et 
al. 2009; Wang et al. 2011). Wang et al. (2007) spec¬ 
ify that halos in overdense environments cannot accrete 
dark matter because of the high velocity dispersion of 
the dark matter. These works generally find that tidal 
effects are strongest for low-mass halos. 

We might expect to see several signatures of this effect 
on our results. First, we would expect slower late-time 
growth of halos in the overdense region compared to ha¬ 
los in the underdense region. Secondly, we would expect 
that within the overdense box, low mass halos would have 
slower late-time growth than high mass halos. Finally, 
we would expect tidal effects to cause the largest differ¬ 
ence between galaxies in the overdense versus underdense 
environment for galaxies with the lowest halo mass. 

In agreement with our first prediction, we see that the 
highest halo mass bin (lO^^'^-lO^^ ® Mq) shows slower 
growth in the Mhaio of the galaxies in overdense regions 
at late t imes (z<l; Figured]). Also, as discussed above 
(Section 14.IL we also find that the formation time of 
galaxies in the higher density environment tends to be 
earlier than in the lower density environment. 

However, when we test the second prediction by com¬ 
paring different mass ranges within the overdense box, 
we do not find the expected trend. As discussed in Sec¬ 
tion |3dIIl the halo age is the oldest for the highest halo 
mass bin in the overdense box, indicating that the late¬ 
time halo growth of low mass halos is larger than that of 
the highest mass halos. If we focus on halo growth from 
z=l to z=0, we find that most mass bins increase by 
about the same factor. A caveat to this result is that the 
smaller-scale environment within the large-scale overden¬ 
sity is not necessarily the same surrounding galaxies of 
different masses. In fact, the number of galaxies within 
2 Mpc of the tracked galaxies decreases with decreasing 
halo mass, so a careful comparison may show that in the 
same smaller-scale density field, lower mass halos grow 
more slowly at late times. 

Finally, against our third expectation, the difference 
between the stellar mass of the lowest mass galaxies 
(lO^^-lO^^'® Mq) is clearly the smallest of any of the 
mass bins, and the similarity of the SMHM ratio sup¬ 
ports that hnding. Higher mass bins (lO^^'^-lO^^'® Mq 
and Mq) have the largest SMHM ratio dif¬ 

ference. 

Therefore, while we find evidence that tidal truncation 
of accretion and/or tidal stripping of dark matter halos 
is occurring in the overdense environment by comparing 
the late-time growth of halos, because the difference in 
formation time is not the largest for low-mass galaxies 
and because the difference in the z=0 SMHM ratio is 
the smallest for low-mass galaxies, we do not think that 
tidal effects are the main driver of our results. 

4.3. Local Galaxy Density 

Tidal effects on galaxies may come from local halos, 
and we have not normalized for the local galaxy density. 
Indeed, in all of the mass bins, the galaxies in the over¬ 
density (G box) have a larger median number of galax¬ 
ies within 2 physical Mpc (pMpc) of the tracked sample 
than the galaxies in the underdensity (V box). In order 
to determine if the local galaxy density within 2 pMpc 
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Fig. 6. — The SMHM ratio (left) and the number of galaxies within 2 pMpc of galaxies (number of nearby neighbors plus the galaxy 
itself), with an additional criterium for the overdense (C box) sample that the z=0 number of nearby galaxies is less than or equal to 2 (3 
including the galax y). Even if we select C box galaxies with lower local densities than that of V box galaxies, the SMHM ratio remains 
higher. See Section l4.3l for discussion. 


is in fact causing the differences in the SMHM ratios, 
we select a subset of C box galaxies in each mass bin 
that ht the additional criteria that at z=0 the number 
of galaxies within 2 pMpc of the C box galaxy is three 
or less, including the galaxy itself, meaning two or fewer 
neighbors. We choose the upper limit of two neighbors 
because the median number of neighbors in the V box 
sample is two in our middle two mass bins ® 

Mq and Mq), one in our lowest mass bin, 

and four in our highest mass bin. Thus we can compare 
this ’’low local-density” C box sample to the complete, 
unchanged, V box sample. 

In Figure [6] we show the SMHM ratio and number of 
galaxies within 2 pMpc in this much reduced sample. 
The small sample size of C box galaxies is most dra¬ 
matic in the highest halo mass bins. In the highest mass 


bin (lO^^'^-IO^^ ® Mq) there is only one galaxy in the 
C box sample until z=1.2, and only two galaxies there¬ 
after. The second highest bin also only has two galaxies 
in the C box (10^^'®-10^^'^ Mq), and the third bin has 
three ® M 0 ). The C box samples in the low¬ 

est mass bin, lO^^-IO^^'® Mq, has 22 galaxies. Due to 
these small samples, the conclusions we can draw are 
necessarily tentative. 

While at early times the C box galaxies always have 
more neighbors, by z=l the numbers are similar in the 
two environments, and by z=0 the number of neighbors 
in the underdense environment is generally slightly more 
than in the overdense environment (by construction, due 
to our selection criteria). We note that this may indi¬ 
cate more mergers in the higher density environment, in 
agreement with Fakouri & Ma (2009). However, even us- 
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ing this specifically selected sample, the SMHM ratio is 
higher in the overdense (C box) sample. Therefore, when 
examining a z<l galaxy sample, simply forcing the num¬ 
ber of galaxies to be the same in 2 Mpc regions is not 
enough to damp out larger-scale environmental effects. 

4.4. Satellites 

We next briefly consider whether satellites could have 
an effect on the SMHM ratio. First, satellites will add 
to the total dark matter mass within one virial radius 
before the stellar particles merge to be identified as a 
single object by HOP, so more satellites might lead to 
lower SMHM ratios. It is also possible that gravitational 
effects from an orbiting satellite could drive gas towards 
the center of the central galaxy and increase the SFR, 
resulting in higher SMHM ratios. 

However, the differences in the number of satellites be¬ 
tween the overdense and underdense samples are quite 
small. None of the central galaxies in the lowest mass 
bin have satellites identified by HOP (which identifies 
stellar groups of 20 particles, or about 2x10^ Mq). In 
the second lowest halo mass bin, the median number of 
satellites is always zero, although the 75 percentile can 
extend to one satellite in both the C and V box sam¬ 
ples. In the second highest mass bin the median number 
of satellites varies between 0 and I, with the overdense 
sample tending to have more. Finally, in the highest mass 
bin, after z~3 the galaxies in the underdense sample have 
a median number of 1-2 satellites while the galaxies in 
the overdense sample have a median number of satellites 
between 0-1.5 (a median of 1.5 can occur with an even 
number of galaxies in the sample). 

While satellites may have something to do with the 
lower SMHM ratio in the underdense galaxies in the high¬ 
est mass bin, it is unlikely that they play a role for the 
majority of galaxies. 

4.5. Gas supply ajfeeted by the large-seale environment? 

We now discuss whether the gas supply available to 
form stars in galaxies of the same halo mass may dif¬ 
fer in different environments. Clearly, the dark matter 
density in the overdense region is higher than in the un¬ 
derdense region. If baryons generally trace dark matter, 
there must be therefore more gas available in the over- 
dense region. Indeed this is the case in our simulation, 
as the global baryon fraction differs by less than 5% in 
the C and V boxes. The question is then: why does more 
gas form stars in halos in the overdense region? 

In recent work, Cen (2014) has found that the number 
of cold streams feeding galaxies correlates with the SFR. 
Therefore, we posit that if galaxies in the overdense envi¬ 
ronment have more cold streams, then more gas may be 
able to form stars. Cen (2014) also argues that streams 
have higher densities in environments that are effected by 
gravitational heating. In our higher density environment, 
individual galaxy halos are more likely to collapse within 
denser filaments or Zeldovich pancakes, which may result 
in more higher density cold streams. 

However, Fakhouri & Ma (2010), examining the dark 
matter-only Millennium Simulation, find that although 
there is more diffuse dark matter around halos in high- 
density environments, the accretion rate of diffuse dark 
matter is higher in low-density environments. They pro¬ 


pose that this may be because diffuse material is dynam¬ 
ically hotter in high-density environments, and therefore 
unable to accrete onto halos. Although determining the 
environmental dependence of galaxy accretion from cold 
streams is beyond the scope of this paper, in future work 
we will compare the cold streams accreting onto galaxies 
in overdense versus underdense regions. 

4.6. Stronger Feedbaek 

We can now discuss the possible effects of including 
AGN feedback or changing the strength of SN feedback 
on our results. 

Including AGN feedback would have the most effect 
on high-mass galaxies (Mhaio >10^^ M©), which are not 
the focus of this work. While there would be some ef¬ 
fect on the galaxies in our highest mass bin, it would 
be quite small (as these halos are <10^^ Mq). Also, 
this effect would be smaller at earlier times, before the 
AGN became massive, although Schaye et al. (2010) find 
that AGN feedback begins affecting the SFR of massive 
galaxies as early as z^i. However, we note that the stel¬ 
lar mass of the C box galaxies in the highest halo mass 
bin is always (to z^6, the entire time that galaxies are 
tracked in the simulation) higher than the stellar mass of 
the V box galaxies. Therefore, we conclude that includ¬ 
ing AGN feedback would not strongly affect our results. 

As discussed above, we posit that much of the differ¬ 
ence in the stellar mass of galaxies in underdense versus 
overdense environments is due to differing gas supply. 
Therefore, decreasing the energy input from supernovae 
feedback would not bring the stellar masses of galaxies in 
these different environments into agreement and might 
increase the difference. On the other hand, increasing 
the strength of SN feedback may bring the stellar masses 
of galaxies in these different environments into better 
agreement. As shown in Figure 12 of Gen (2011c), the 
specific cold gas inflow rate is always higher than the 
sSFR, in both the C and V boxes, which is likely due 
to SN feedback regulating star formation rates. There¬ 
fore, it may be possible to eliminate any differences with 
environment if we were to increase the strength of our 
feedback to lower the SFR at all halo masses. As in¬ 
creasing the feedback strength would affect halos in both 
the overdense and underdense environments, we expect 
that a small change would not be enough to bring the 
stellar masses into agreement. Therefore we are hesitant 
to invoke increased feedback to remedy this difference, 
as in several papers we have compared our simulations 
to observations and found good agreement. For exam¬ 
ple, the SFRD in our two extreme regions brackets the 
observed global SFRD (Gen 2011c), and dramatically af¬ 
fecting the SFR of galaxies in our simulations would have 
a significant effect on this measurement. 

As shown in Cen (2011c), the specific cold gas inflow 
rate is higher than the sSFR for galaxies in both the C 
and V boxes, indicating self-regulation of SF. However, 
the difference between the specific inflow rate and the 
sSFR is highest at high redshift (z>3), and at these high 
redshifts V galaxies have a higher inflow to sSFR ratio. 
Therefore, we find that SN feedback is making SF less 
efficient in galaxies in the underdense environment. Be¬ 
cause our current feedback scheme seems to affect the gas 
consumption rate in V galaxies more strongly than in C 
galaxies, we infer that increasing the feedback strength 
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would only exacerbate the difference in the populations. 

Although we cannot use SAMs to make strong predic¬ 
tions for hydrodynamical simulations, we note that in 
Jung et al. (2014) the authors vary their star formation 
feedback strength by a factor of four and find no qual¬ 
itative difference in their result for low mass halos that 
galaxies in dense environments have higher stellar masses 
than galaxies in underdense environments. 

Finally, even if we devised a feedback prescription that 
eliminated the differences in SFR and therefore stellar 
mass in halos in different large-scale environments, our 
results are interesting because they indicate that gas ac¬ 
cretion differs across different large-scale environments, 
and is not solely dependent on halo mass. 

4.7. How environment affeets the SMHM ratio 

We conclude that the SMHM ratio depends on more 
than the halo mass or local (2 pMpc or less) galaxy den¬ 
sity at z=0. Our results point to a picture in which the 
combination of earlier formation times, enhanced inter¬ 
actions, and more filaments are the likely causes of the 
increased SMHM ratios in galaxies in higher density en¬ 
vironments. 

First, formation time plays a role. Galaxies in the over- 
dense region have earlier formation times. An earlier for¬ 
mation time means that the galaxies in the overdense re¬ 
gion grow more quickly at earlier times, when all galaxies 
have higher SFRs and more star-forming gas available. 
Therefore, their higher halo mass at z>l allows them to 
form more stars when rates are high. As we are con¬ 
sidering halos that all have z=0 masses below 10^^ Mq, 
the halos are not massive and hot enough to quench cold 
accretion and/or star formation at early times. 

Second, interactions with nearby neighbors play a role. 
Galaxies in overdense environments have more nearby 
neighbors, particularly at early times. In Figure [6l where 
we have forced the z=Q number of neighbors in the G box 
to be less than in the V box, we see that at early times 
the SMHM ratio is lower and the number of neighbors 
is higher in the galaxies in the overdense region than in 
those in the underdense region. The SMHM ratio in¬ 
creases in the C box galaxies relative to the V box galax¬ 
ies, as the number of neighbors decreases in the C box 
to equal the number of neighbors are V box galaxies. 
As discussed in Gen & Safarzadeh (2014), interactions 
first drive enhanced star formation in massive galaxies, 
then at later times affect lower mass galaxies. In agree¬ 
ment with this scenario, in Figure [6] (and Figure [3]) the 
SMHM ratio of the C box galaxies crosses that of the V 
box galaxies at later times for lower mass galaxies. 

Finally, filamentary structure affects galaxy growth. 
As we discuss in Section l45l Gen (2014) finds that num¬ 
ber of filaments correlates with SFR. Overdense regions 
are likely to have more filaments that can feed galaxies 
through smooth accretion and increase SFRs. In par¬ 
ticular, the growth in both the stellar and dark matter 
mass of galaxies in the overdense region happens at ear¬ 
lier times, when filamentary accretion is more likely to 
be able to strongly increase the mass of galaxies, and gas 
consumption into stars is more efficient (Cen 2011c). At 
later times, filamentary accretion becomes less important 
in general, which is why filaments can impact the SMHM 
ratio even of galaxies in the large scale overdensity that 
are in local underdensities by z=0. 


The relative SMHM ratios of galaxies in different en¬ 
vironments depends on the mass of the galaxies. In the 
two most massive bins in Figure|3]the difference between 
the SMHM ratio in the two environments peaks before 
z=0. This is because gas accretion and star formation 
have slowed in the high mass galaxies in the overdense 
environment, but are still continuing in the high mass 
galaxies in the underdense environment. The difference 
in the SMHM ratio remains similar from z=l to z=Q in 
the second mass bin (10^^ ®-10^^ ® Mq), and is increasing 
for galaxies in the lowest mass bin. We see this trend with 
mass reflected in Figure [3 that the sSFR in the highest 
mass bin is lower in the overdense than in the underdense 
environment, while at the lowest halo masses the sSFR 
is similar in the overdense environment. Therefore, we 
suspect that for even higher masses, by z=0 the differ¬ 
ence in SMHM ratio would be entirely wiped out by the 
higher sSFR in the low density environment. 

In summary, the environment affects the formation 
time of halos, which in turn affects how gas-rich the uni¬ 
verse is when they are massive enough to quickly accrete 
gas and form stars. Also, galaxies in higher density large- 
scale environments have a larger number of neighbors, 
particularly at early times. These interactions can drive 
gas instabilities and star formation. Finally, galaxies in 
higher density environments are likely fed by more, and 
higher density, filaments. 

5. COMPARISON WITH OTHER WORK 

Much of the work performed examining how environ¬ 
ment, often measured through clustering, affects galaxy 
properties has been performed using SAMs. A dark 
matter-only N-body simulation is used and galaxies are 
assumed to populate halos from early times. The assem¬ 
bly history of the halos, such as smooth accretion and 
mergers, is also assigned to the galaxy. For example, 
the accretion of gas follows the accretion of dark mat¬ 
ter using the global baryonic fraction, then cools at a 
prescribed rate until it can form stars. The internal pro¬ 
cesses of galaxies, such as gas cooling, star formation and 
feedback processes, are generally simplified prescriptions 
that are tuned so that the z=0 population matches ob¬ 
served galaxy populations (e.g. as described in Jung, Lee 
& Yi 2014). 

Using such a SAM, Jung et al. (2014) find a small dif¬ 
ference in regions of different large-scale (7 h“^ Mpc) 
density, specifically that low-mass (<10^^ Mq) halos 
have slightly higher stellar masses in high-density envi¬ 
ronments. 

Our results qualitatively agree with those of Jung et 
al. (2014), although they find a much smaller difference 
in the stellar mass of galaxies in overdense versus under¬ 
dense regions, and only for low-mass halos (lO^^-lO^^'^ 
Mq/H). They show that in their simulation this is be¬ 
cause there is more cold gas in low-mass halos in higher 
density environments. We have also posed this as an 
explanation for our different stellar masses. However, 
we see a much larger effect for a much broader range 
of galaxy masses. We consider a few possible explana¬ 
tions. First, Jung et al. (2014) considers environment 
on the scale of 7 h~^ Mpc while we are considering en¬ 
vironment measured on ^25 Mpc scale. Although Jung 
et al. (2014) tested that scales from 3-9 h~^ Mpc did 
not effect their results, Gao et al. (2005) found that as- 
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Fig. 7. — The halo sSFR (SFR/M£)jvfh,aio) ^ function of dark matter halo mass. Each point is the median value of 50 galaxies binned 
by halo mass. The shaded region shows the 25-75% range of values. The red denotes the C box galaxies, and the blue denotes the V box 
galaxies. We overplot results from Crain et al. 2009 in the z=0 panel for comparison. From left to right the redshifts are 2=3.1, 1.0, 0. 
The halo sSFR tends to be higher for galaxies in the overdensity (C box) than in the underdensity (V box). As the redshift decreases, the 
difference between C box galaxies and V box galaxies become more dramatic. See discussion in Section [5] 


sembly bias increases with increasing length scale, and 
Jung et al. (2014) do not extend to the scale we use to 
define environment. Secondly, we directly calculate the 
cooling rate using the gas chemistry, and perhaps gas 
is more metal-rich in the overdense environment and so 
can cool more quickly. Indeed, Cen (2013) finds that gas 
flowing into halos (with negative radial velocities) tends 
to have higher metallicity in the overdense C box than 
in the underdense V box. It is possible that metals can 
escape lower density filaments more easily therefore re¬ 
sulting in less metal-enrichment of galaxies in low density 
environments. This is an extension of the discussion of 
metal enrichment of clusters via filaments by Rasmussen 
& Ponman (2009). However, we suspect that the most 
important difference is that we directly calculate any hy- 
drodynamical effects from neighbors. As we discuss in 
Section l4?7l we believe that one driver of our results is 
interactions with nearby galaxies (most of which are not 
satellites) driving gas instabilities and inflows that in¬ 
crease SFRs. This is a process that is not treated in 
SAMs. 

We can also compare with the other hydrodynamical 
simulation that compares galaxies in different environ¬ 
ments, and find that our results differ from those in Crain 
et al. (2009). The GIMIC simulations use GADGETS 
to hydrodynamically “re”-simulate regions of the Millen¬ 
nium simulation at (-2, -1, 0, -1-1, -l-2)tT (on the scale 
of ^20 Mpc) of the mean density to determine what 
drives differences in the SFRD. While the peaks of their 
curves are between z=2-3, in agreement with our results, 
they find a smaller difference in the normalization of the 
SFRD curves in the high- and low-density environments. 

Other properties of their galaxy populations differ from 
ours. For example. Grain et al. (2009) do not find much 
difference in either the shape or the magnitude of the 
stellar mass functions in their over dense versus under 
dense regions. However, while they reproduce the over¬ 
all number of galaxies more massive than 10® Mq, their 
stellar mass functions do not match the observed shape 
of the z=0 stellar mass function of Li & White (2009) (as 
discussed in their Figure 3). We note that this group has 
since run the ‘Evolution and Assembly of GaLaxies and 


their Environments’ (EAGLE) simulation project, which 
has an excellent fit to the z=0 galaxy stellar mass func¬ 
tion (Schaye et al. 2015). However, the EAGLE simula¬ 
tions do not compare different large scale environments, 
so we cannot compare our work to these simulations. In 
Figure 3 of Cen (2011) the luminosity function of galaxies 
in the combined C and V boxes is shown to be a good fit 
to Blanton et al. (2003), although the highest luminos¬ 
ity, i.e. most massive, galaxies require a post-processing 
addition of AGN feedback to match observations. In ad¬ 
dition, Kreckel et al. (2011) focus only the V box, and 
reproduce the luminosity function of void galaxies found 
by Hoyle et al. (2005). The void luminosity function dif¬ 
fers from the total luminosity function both in that the 
normalization of the low luminosity end drops by more 
than an order of magnitude and the knee of the function 
shifts. 

In this paper we track galaxies over time rather than 
look at the galaxy population at individual redshifts. 
However, in order to determine whether we would get 
the same results, in Figured we make a halo sSFR fig¬ 
ure similar to Figure 8 in Grain et al. (2009). Here we 
look at 2=3.1,1.0, and 0. We bin our galaxies by halo 
mass, with 50 galaxies in a bin, except the highest mass 
bin has between 50-100 galaxies. In red we plot the G 
box galaxies and in blue we plot the V box galaxies. 
The solid lines and symbols are the median values and 
the shaded region denotes the central 25-75% of values. 
We briefly note that the low mass galaxies at 2=0 have 
masses below our lowest mass bin in Figures [T]l6l and 
most of these galaxies are not tracked even to 2=1. The 
precipitous drop in the sSFR is quite probably due to 
poor resolution of these galaxies, or may possibly indi¬ 
cate that the environment is dramatically affecting the 
SFRs of these galaxies beyond two virial radii of their 
nearest neighbor or that these are splashback galaxies. 
The dashed lines at 2=0 show some of the Grain et al. 
(2009) data. First, our halo sSFR agrees well with the 
high mass end of the halo sSFR in Crain et al. (2009), 
but is relatively flat with halo mass and does not drop 
at lower masses. We speculate that this is because we do 
not include kinetic feedback as in the GIMIC simulations. 
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so we do not have a sudden drop below about 10^^ Mq 
from gas blow-out. Unlike Crain et al. (2009), we find 
that at lower redshifts the halo sSFR is clearly different 
in differing environments, across all of the halo masses 
we probe. Even at 2=3.1, we see a hint that the halo 
sSFRs of galaxies in the overdense environment is larger 
than that of galaxies in the underdense environment. At 
earlier redshifts, the sSFRS in the two environments are 
indistinguishable. 

The reason for these differences is difficult to determine 
as the simulations are quite different. Crain et al. (2009) 
use GADGET3, a smooth-particle hydrodynamics code. 
Their gravitational softening length is fixed in physical 
space at 2<3 to 1 h~^ kpc in the intermediate resolution 
runs (from which most of their results are drawn), which 
is twice our resolution scale. Their cosmology is matched 
to the Millennium Simulation, so is also different from 
ours. As discussed above, they also employ kinetic rather 
than thermal feedback. 


6. CONCLUSIONS 

In this paper we have used a cosmological hydrody- 
namical simulation to examine the stellar mass to halo 
mass ratio of central galaxies in two different large-scale 
(^20 Mpc) environments. When we focus on the growth 
of dark matter halos (Figure [T]), we reproduce the as¬ 
sembly bias found by Gao et al. (2005)-galaxy halos 
in overdense environments have earlier formation times 
than those in underdense environments. 

Importantly, we have found that the halo mass alone 
cannot determine the SMHM ratio of galaxies. Specifi¬ 
cally, in the large-scale overdensity central galaxies with 
halo masses from ® Mq have higher SMHM ra¬ 

tios than galaxies in a large-scale underdensity (Figure 
El). Even when we force the number of galaxies within 2 
pMpc at 2=0 in the overdense regions to be similar to or 
smaller than the number of surrounding galaxies in the 
underdense region, we cannot erase this difference. 

We posit that this result is due to the earlier forma¬ 
tion times of halos in a large-scale overdensity and to the 
larger number of nearby galaxies and filaments at early 
times. Because the halos in overdense environments are 
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more massive at early times, they can accrete gas and 
therefore form stars more quickly than those in under¬ 
dense environments. Although at later times the V box 
halos grow more quickly than the C box halos, by that 
time there is generally less gas accretion, so SFRs are nec¬ 
essarily lower and the SMHM ratio does not equalize. In 
addition, particularly at early times, halos in overdense 
environments have a larger number of neighbors, which 
can drive instabilities in the gas and increase the SFR of 
central galaxies. Einally, there are more filaments feed¬ 
ing galaxies in higher density environments, particularly 
at early times. 

We find that the 2=0 sSFR of galaxies is higher in 
the underdense region, in agreement with observational 
results (Rojas et al. 2004; Grogin & Geller 1999,2000; 
Szomoru et al. 1996). To test our results more ro¬ 
bustly, we need more detailed observations of the SFR 
and stellar masses of a large sample of galaxies in deep, 
large voids. In particular, more comparisons of the Tully- 
Fisher relation using red bands to determine stellar mass 
would be useful. One difficulty will be to make sure we 
are comparing only central galaxies in both underdense 
and overdense environments. 

Different SMHM ratios in different large-scale environ¬ 
ments means that assigning stellar masses based entirely 
on halo masses in simulations may not produce a realistic 
galaxy population, and including the 2=0 local density in 
the stellar mass assignment will not compensate for the 
different environments in which the galaxies previously 
evolved. 
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